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Abstract] The Hessian fly Mayetiola destructor Say and the barley stem gall midge Mayetiola hordei Kieffer are two 
sibling phytophageous species causing[] annually[] important losses to cereals in Tunisia. Damage to wheat has usually 
been attributed to Mayetiola destructor but there is no strict association between the cereal specied] wheat or barley[] 
and the Mayetiola species[] destructor or hordei[]. An effective pest management program requires an accurate 
analysis of the midges’ genotypes. In the present study[] we used the RAPD technique in conjunction with mating 
bioassays and mitochondrial DNA typing[] to assess the extent of genetic variability and taxonomic relatedness[] within 
a Mayetiola population infesting wheat in the North of Tunisia. Phylogenetic analysis inferred from the RAPD showed 
that the studied population had a variable genetic background that could be due to a complex epidemiological situation 
where wheat samples would be co-infested by both M. destructor and M. hordei. Although mating assays did not 
result in oviposition in a few cased] 2/14[[] mitochondrial typing based on the Cyt b gene restriction analysis showed 
that all samples were M. destructor ones. Our results are additional evidence that RAPD variability is not a reliable 
criterion for taxonomical inference and provide support for allopatric distribution of M. destructor and M. hordei 


species[] in Tunisia. 
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Dipterans of the Mayetiola genus are destructive 
pests that cause important losses on cereals. The 
insect' s female lays 200 to 300 eggs on the leaves of 
infested plantd] Ratcliffe and Hatchett 1997[]. Upon 
hatching[] larvae enter the leaf sheath and feed near the 
crown or at nodes[] resulting in stunting of the plant and 
development of a dark green colour in the leaves 

[] Beyers and Gallun[] 1972[] Zantoko and Schukle[] 
1999 []. The genus 
phytophageous species[] namely Mayetiola destructor 


Mayetiola includes two 
Say[] commonly known as the Hessian fly[] and 
Mayetiola hordei Kieffer[] commonly known as the 
barley stem gall midegel ]. 

There is evidence that Mayetiola destructor 
historically followed wheat from its original habitat in 
South-West Asia until it dispersed throughout most 
wheat-growing areas in the world[] Briggle et al .[] 
1982[]. Nowadays[] this species is found in Asia[] 
Siberia[] Europe[] North Africa[] North America and New 
Zealand] Ratcliffe and Hatchett[] 1997[]. Wheat is the 
preferred host of the Hessian fly[] although barley] ryel[]] 


triticale and several volunteer grasses [] such as 





Agropyron and Aegilops[] may occasionally serve as 
hosts when wheat is not availabl¢] Buntin and Raymer[] 
1989[] Ratcliffe and Hatchett[] 1997[]. 

Mayetiola hordei was mentioned only in North 
Africa and Southern Europd] Skuhrava et al .[] 19840 
Gagne et al.[] 1991[]. It reproduces on barley but it 
can also sometimes be found on wheat[] Mezghani et 
al.[] 2002a[]. Thus[] there is no strict association 
between the cereal species] wheat or barley[] and the 
Mayetiola speciei] destructor or hordei]. 

The two phytophageous species of the Mayetiola 
genus can difficultly be distinguished morphologically . 
In fact[] Gagne et al [] 19911] using scanning electron 
microscopy to analyse midges from Morocco[] reported 
differences in the shape of female postabdomen[] the 
structure of male terminalia and the number of spicules 
covering the puparium. In Tunisia] Makni et al. 

O 2000[] identified diagnostic allozyme markers of 

taxonomic value for the two species[] however[] this 
technique requires refined conservation conditions to 
preserve allozyme activity. 


More recently[] Mezghani et al. [] 2002b[] 
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investigated the correlation between allozyme markers 
O Makni et al.[] 2000[] and mitochondrial DNA 
sequence variation[] by sequencing the cytochrome b 
O Cyt b[] gene from a collection of Tunisian M. 
destructor and M. hordei samples. The authors 
reported three mitochondrial haplotypes namely 
haplotype A found in M. destructor and haplotypes B 
and C found in M. hordei. Based on this sequencing 
data[] it was possible to diagnose the two sibling species 
of Mayetiola by means of the Restriction Fragments 
Length Polymorphisn{] RFLP[] of the Cyt b gend] 800 
bp[]. In fact[] using the VspI endonuclease[] Mezghani 
et al.[] 2002c[] showed that haplotype A[] in M. 
destructor[] was represented by two restriction fragments 
[] 540 and 260 bp[I[] whereas haplotypes B and C[] in 
M. hordei[] were represented] respectively[] by three 
[] 540[] 170 and 90 bp[] and two] 630 and 170 bp[] 
restriction fragments. Although this is a simple and 
reliable approach[] inferences regarding taxonomical 
specificity need to be further consolidated by conducting 
mating bioassays[] according to the biological species 
concept] Mayr{] 1974[]. 

In Tunisia[] widespread outbreaks of Mayetiola[] 
causing extensive crop losses on cereals[] have been 
reported since the early 1900' s[] Pagliano[] 1935[] 
Miller et al .[] 1989[] Makni[] 1993[]. Particularly on 
wheat[] damage caused by these insects was estimated at 
40% of durum wheat[] Triticum durum  Desf.[] and 
26% of bread wheal] Triticum aestivum L.[[] Makni[] 
1993[]. Breeding Tunisian wheat varieties for genetic 
resistance to Mayetiola has been a major research focus 
in order to reduce yield losses[] Bouktila et al.[] 
2005[]. However] for this approach to succeed[] 
Mayetiola populations occurring on wheat should be 
studied in terms of genetic variability and species 
identity. In this context[] the Randomly Amplified 
Polymorphic DNA[] RAPD[] technique is a rapid and 
cost-effective methodology that has successfully been 
used to detect polymorphism in several insect species 

O Black et al.[] 1992[] Vanlerberghe-Masutti and 
Chavigny[] 1998[] Naber et al .[] 2000[]. 

Our objectives in the present study were[[] i[] to 
assess the extent of genetic variability within a 
Mayetiola population infesting wheat in one of the 
Tunisian cereal-growing regions and[] ii[] to clarify 
whether the genetic community of Mayetiola on wheat[] 
in the studied region[] is con-specifid] all belonging to 
M. destructor[] or made by a mixture of M. destructor 
and M. hordei . 


1 MATERIALS AND METHODS 
1.1 Insect material and rearing methods 


A hundred wheat plants[] infested with Mayetiola 
pupad] were randomly sampled in April 2004[] from 


private fields located along a 5 kilometers Northwest / 
Southeast axis[] across the region of Jedeida[] 30 km to 
the Northwest of Tunid] Fig. 1[]. Those infested plants 
were transferred to the laboratory[] inside plastic bags. 
The rearing method was similar to that described in 
Cartwright and Lahud] 1944[[] Mayetiola larvae at the 
puparium stage were conserved inside cheesecloth- 
covered cages in a greenhouse where conditions are 
favourable for adult emergence] 12 h photoperiod from 
6[ DO to 18[00[] 20 + 1°C temperature and 60% relative 
humidity[]. On this regime[] adults emerged from 5[B0 
to 8[]30] McKay and Hatchett[] 1984]. Thus[] virgin 
adults were obtained readily by inspecting the cages at 
regular 2 min intervals from 5[]30 to 8030 and 
determining the sex. Eighteen adult insects having 
emerged in these conditions were used for RAPD 
analysis[] whereas 28 virgin others[] 14 males and 14 
females[] were used for mating tests and subsequent 


PCR-RFLP analysis. 
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Fig. 1 A map showing the geographic location of the 


studied population of Mayetiold] open circle[] and 
the sampling axid] arrow[] 


1.2 DNA isolation 
For molecular analyses[] we isolated DNA from 
each insect individual following the Doyle and Doyle 
[] 1987[] DNA extraction protocol. Each fly was placed 
in a separate Eppendorf tube and ground with a 
miniature pestle in 200 pL of extraction buffet] 100 mM 
Tris-HCl pH 8[] 20 mM EDTA[] 1.4 M NaCI[] 2% 
CTAB and 0.2% 2-mercaptoethanol[]. After incubation 
at 65 ^C for 30 min[] 200 pL of chloroform[] isoamyl 
alcoho[] 24: 1[] were added[] mixed and centrifuged at 
8000 rpm for 5 min at 4 °C and the supernatant was 
transferred to a new tube. Nucleic acids were 
precipitated overnight with isopropanol] 1:1[] at - 20 
°C. The mixture was centrifuged at 13500 rpm for 15 
min at 4 C[] the supernatant discarded and the 
remaining pellet washed in 70% Ethanol for 10 min and 
dissolved in 20 uL sterile water. 
1.3 RAPD amplification 
DNAs isolated from 18 Mayetiola individuals were 
amplified in a volume of 25 „LO Williams et al .[] 
1990[] with 25 ng of template genomic DNA[] 1.75 units 
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of Taq DNA Polymerasd] QBiogene[] France[ T] 100 yM 
of dNTPs[] 100 ng of primer and 2.5 uL of 10X reaction 
buffer supplied by the enzyme manufacturer. A set of 5 
primers of arbitrary sequences[] supplied by Operon 
Technologied] USA[T] was used for RAPD amplification 

[] Table 1[]. Amplifications were performed in a thermal 
cycler’ Crocodile [I| '[] Qbiogene[] France[] programmed 
as follows[] 1 cycle of 5 min at 95 ?C[] 35 cycles of 30 
seconds at 95 °C[] 1 min at 35 °C and 1 min at 72 C[] 
and a final cycle of 5 min at 72 °C. A negative control 
without DNA was added in each run to test for 
contamination and each reaction was replicated twice to 
check for reproducibility of the RAPD patterns. 


Table 1 Codes and sequences of primers used[] total number 
of bands analysed and fragment size 





Total Band 

Pri S number size 

Ng i alae of bands range 

generated O bpl] 
OPD08 5'-GTGACGTAGG-3' 22 250 - 1 800 
OPD10 5'-GTGATCGCAG-3' 22 300 — 3 000 
OPDI5 5'-TTCCGAACCC-3" 19 220 - 2 800 
OPDI6 5'-AGCCAGCGAA-3' 30 220 - 2 500 
OPD20 5'-GTTGCGATCC-3' 19 220 - 1 800 
Total - 112 220 - 3 000 

Mean/primer = 22.4 - 
Amplification products were separated by 


electrophoresis in a 1.5% agarose gel in 0.5% TBE 
buffer. Amplification products were visualized under 
UV lights and photographed after staining the gel in a 
solution containing 1 pg/mL ethidium bromide. The 
molecular weight of each DNA band was estimated by 
comparing with a co-migrating 1 kb ladder. 
1.4 RAPD data statistical analysis 
Under the assumption that variation in banding 
patterns represents allelic segregation at independent 
loci] each RAPD band was treated as a two-allele 
system coded as l[] presence of the band[] or O 
[] absence of the band[]. The extent of genetic variability 
within the studied population was estimated by 2 
parameterd] [] i[] the percentage of polymorphic 
markers] and] ii[] the genetic distancd] Nei and Lil] 
1979[] calculated from the binary data for all pair-wise 
combinations of samples using the formulal] D = 2 
nap] 2 nag + n4 + ng[T [] where ngis the number of 
bands present in individual A only[] ng the number of 
bands present in individual B only and nag the number 
of bands present in both A and B. Using this Nei and 
Li' d] 1979[] distance[] cluster analysis was performed 
to study genetic relationships between Mayetiola 
individuals. Two methods of cluster analysis[] the 
neighbour-joining]] NJC] Saitou and Neill 1987[] and 
the unweighted paired-group method for the arithmetic 
averagd] UPGMA[]] Sneath and Sokal[] 1973[[] were 
used. The" RAPDistance Programs Version 2.00” 


[] Armstrong et al.[] 1999[] was used for genetic 
distance matrix calculation and NJ dendrogram 
construction[] whereas the" PHYLIP package Version 
3.572 C"[] Felsenstein] 1995[] was used for the 
UPGMA dendrogram construction . 

1.5 Mating bioassays 

The mating assays were carried out in a 
greenhouse set at 12 h photoperiod from 6[00 to 180000 
20+1 °C and 60% relative humidity. Five seeds of a 
Mayetiola-susceptible wheat cultivar] " BlueBoy " 

[] kindly provided by the Kansas State University[] 
Manhattan[] USA[T[] were planted in each of 14 soil- 
filled pots. In each pot[] 1 virgin male and 1 virgin 
female of Mayetiola[] having recently emerged from 
infested wheat samples collected from Jedeida[] were 
released for mating inside a cheesecloth tent placed over 
the pot. After 48 h[] tents were removed and the 28 
dead adult flies were collected. Plants were examined 
individually under a microscope for eggs. 

1.6 Mitochondrial DNA amplification 

The 28 Mayetiola individuals used for mating tests 
were subjected to mitochondrial DNA amplification. An 
800-bp DNA segment of the cytochrome H] Cyt b[] gene 
was amplified using a couple of primers CPI[] 5'- 
GATGATGAAATITTGGATC-3[] and CP2 [] 5*- 
CTAATGCAATAACTCCTC-3'[][] Mezghani et al.[] 
2002c[]. Each PCR reaction mixture contained 10 ng of 
template DNA[] 2.5 pL of 10X reaction buffer 

[] Qbiogene[] France[T] 1 unit of Taq DNA Polymerase 

O Qbiogene[] France[ T] 100 mM of dNTPs[] 0.7 uM of 
each primer. The final reaction volume was 25 pL. 
Thermal-cycling was performed in a" Crocodile [I| " 
PCR machind] Qbiogene[] France[T] as follows[] a cycle 
of 1 min at 95 *C[]35 cycles of 1 min at 95 ^C[] 1 min 
at 48 ^C and 1 min at 72 ^C[] and a final cycle of 5 min 
at 72 C. 

1.7 Restriction analysis 

15 pL of each PCR product were separately 
digested] overnighti[] by 2.5 U of the restriction 
endonuclease Vsp [J] ATI[]JAAT[] at 37 *C[] Mezghani 
et al .[] 2002c[]. The digestion products were separated 
by electrophoresis in a 296 agarose gel stained with 
ethidium bromide and visualised under UV lights. 


2 RESULTS 


2.1 RAPD analysis 
112 different RAPD bands were generated by the 5 
primers analysed. The total number of bands scored per 
primer ranged from 19[] OPD15 and OPD2O[] to 30 
[] OPDIG[] with an average of 22.4 bands per primer. 
The size of the amplified fragments ranged from 220 to 
3 000 base pair$] bp[]. Table 1 summarises these data. 
Many different patterns were detected within the 
18 studied Mayetiola individuals from Jedeida 
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population[] indicating heterogeneity within the 
population. Fig. 2 shows an example of RAPD profile 
ilustrating this variability. Among the 112 bands 
scored] 101 ones [] 90.196[] were polymorphic. 
Furthermore[] the studied individuals showed differences 
in the levels of genetic distance[] ranging from 0.06 to 
0.63[] with a mean genetic distance equalling 0.27 
O Fig. 3[]. This indicated that the studied samples[] 
within the Tunisian population of Mayetiola infesting 
wheat[] had a variable genetic background. 

Cluster analysis based on the genetic distances 
yielded two similar trees by both the UPGMA and the 
NJ methods[] Fig. 4[]. Both analyses grouped the 
population into 2 distinct clusters[] a major cluster with 
14] NJ[] or 16[] UPGMA[] individuals and a minor 
cluster with 3] NJO or 2] UPGMA[] individuals. The 
minor cluster individuals were outside the remaining 
Mayetiola samples reflecting an evolutionary 
tendency[] possibly with respect to their taxonomical 
specificity . 


10 11 12 19 14 15 16 17 18 





Fig. 2 RAPD pattern generated by using DNA from 
18 Mayetiola individuals of the Jedeida population[] 
with primer OPDIO 
LJ Molecular size marker{] C[] Negative control without DNA. 
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Fig. 3  Pair-wise genetic distanced] Nei and Li[] 1979[] 
between 18 analysed Mayetiola individuals infesting 
wheat in Jedeida 


2.2 Mating assays and cytochrome b gene 
polymorphism 

Plants examination[] 48 h after infestation[] showed 
that wheat plants of the" BlueBoy" cultivar harboured 
Mayetiola eggs in 12 different assays[] indicating a 
fecundity rate of 85.796[] 12/14[] within the studied 
Mayetiola samples having emerged from infested wheat 
in Jedeida. 


The 14 pairs of Mayetiola used in the mating 
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Fig. 4 Cluster analyses based on Nei and Li' 4] 1979[] 
genetic distance[] showing genetic relationships between the 
18 individuals of Mayetiola collected from wheat in Jedeida 

AQ] Neighbour-Joining dendrogram[] BL] UPGMA dendrogram. 









































































































































study were subjected to DNA isolation followed by 
amplification of the Cyt b gene and its restriction by 
Vsp I. This digestion yielded a uniform pattern 
consisting of two fragments sized 540 bp and 260 bp 
O Fig. 5[[] indicating that all the analysed individuals 
are M. destructor oned] haplotype AL. 
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Fig. 5 RFLP pattern of PCR-amplified cytochrome H] Cyt bO 
genel] digested by the restriction endonuclease Vspl 
LJ] Molecular size marker. Lanes 1 — 28[] Digestion of Cyt b 
amplicons from the 28 Mayetiola individuals 
used in mating tests shows a uniform pattern 
corresponding to haplotype A of 
Mayetiola destructor . 


3 DISCUSSION 


The aim of this work was to study a Mayetiola 
population occurring on wheat in Jedeida[] North 
Tunisia] in terms of genetic variability and species 
identity. Using the RAPD technique[] we showed a 
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substantial amount of variation in RAPD profiles among 
Mayetiola individuals[] indicating a variable genetic 
background of insects from this genus that infest wheat 
plants. The dendrograms generated from the genetic 
distance matrix by the UPCMA and NJ methods showed 
only minor differences[] between both methods[] in the 
grouping of the Mayetiola samples thus it was 
concluded that the analysed population derived from two 
founder groups[] founder group 1 comprised the majority 
of analysed flies[] whereas founder group 2 included a 
few number[] 3 by the NJ method and 2 by the 
UPGMAL]. Previous studies have already shown that the 
RAPD technique can be an indictor for species 
identificatiof] Guizani et al .[] 1994[]. Thus[] the high 
genetic heterogeneity within Mayetiola individuals 
infesting wheat[] coupled with the uniqueness of a minor 
group of individuals in both cluster analyses could 
correspond to the existence of two sibling species of 
Mayetiola infesting wheat in Tunisia] namely M. 
destructor[] having a natural preference for wheat[] and 
M. hordei] which could occasionally reproduce on this 
plant] Mezghani et al .[] 2002a[]. 

In order to clarify whether the genetic community 
of Mayetiola on wheat[] in the studied region[] is 
conspecifid] all belonging to M. destructor[] or made 
by a mixture of M. destructor and M. hordei[] we 
performed mating bioassays between 14 couples of 
Mayetiola having emerged from infested wheat in the 
same region. Although two out of the fourteen achieved 
assays did not result in oviposition[] the mitochondrial 
typing of crossed insects indicated that they all belonged 
to Mayetiola destructor species. 

Although the individuals were different from each 
other with respect to randomly amplified fragments 
profiles[] they all belonged to a unique species[] as 
inferred from the Cyt b gene restriction analysis. This 
comparison indicated that the RAPD profiles were not 
associated with the species trait[] in the analysed 
population of Mayetiola . Indeed[] inferences regarding 
taxonomical specificity and/or relatedness could not be 
ascertained by the RAPD] Guizani et al .[] 2002[[] as 
the mechanisms underlying product amplification remain 
obscure[] due to the randomly nature of this technique 

O Williams et al .[] 1993[] Mori et al .[] 1999[]. Thus[] 
the high extent of genetic diversity within the analysed 
Mayetiola destructor samples] estimated by the 
percentage of polymorphic bands and the Nei and Li’ s 

[] 1979[] pair-wise genetic distances[] would rather be 
attributed to[[] if] the ancient colonization of Hessian fly 
on wheat in Tunisia[] in a co-evolution process of the 
pest with its preferred host plant. In fact[] Mayetiola 
destructor would have followed wheat from it’ s original 
habitat in South West Asia and dispersed to North 
Africd] Briggle et al .{] 1982[] Naber et al .[] 200000 


and] iil] the sexual reproduction regime that makes the 


population both heterozygous and  heterogeneous[] 

contrarily to some insect species[] such as aphids[] were 

populations tend to constitute parthenogenetic clones 
[] Puterka et al .[] 1993[]. 

The absence of eggs in two bioassays does not 
imply the existence of a mating barrier but is possibly 
the consequence of a non-preference for oviposition[] 
conditioned by the plants vigour in the corresponding 
pots. Indeed[] the plant genotype and vigour are 
essential factors that influence oviposition preference in 
many phytophageous insect species[] Price[] 1994[] 
Larsson et al .[] 1995[]. 

The results obtained in the present study raise the 
hypothesis that the two species of M. destructor and 
M. hordei would be separated geographically 

[] allopatric[] in Tunisia[] with the presence of only M. 
destructor in Northern regions such as Jedeida[] in 
association with extensive and mono-variety wheat 
culture. This suggests that Mayetiola-controlling efforts 
should be directed towards the resistance to M. 
destructor in the North of Tunisia and M . hordei in the 
SoutH] Makni[] 1993[]. A more extensive sampling will 
help study the insect populations occurring in Tunisia 
and analyse their genetic relationships. 
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unmiu HW DLL. Mayetiola destructor Say L] UUUUUU M. hordei Kieffer UL D] E] BL D] D. EE B. W. DE. EE HE. HDI DE Ul 
uüuggusubBaübBggsaiggsusildgglsilil d d DOO CELO CD DL BL O CELO C] OEC] ELO ELO DOE U D EE C 
HU D OU destructor [] hordei] DU D UU D E. B D. D B HL D. DEB] D. D. EE ODE HE DI DE DI LI uuuiWilitlhibitlLibLl 
guuggggdaggudiyuguulüult DNA RAPo Hi HH] D] 0 H H] D HH C] I pPNAD FH GO D] CD. HU 7T. D. E DIU. DI I 
HUHUHUUUUUUUUUUUUUUUUUUUUUUUUUUUUU RAPDDCD CLO DO DO D U D CDI 
HUHHOUUUUUUUUUUUUWMUUUUUUUUUUUUUUUUUUUUUUUUUUuUUUid 
uguuguguuuuudiagddteismumnaunmuauagaggu»bsapaaaggagaputdsadtgbtüdedgctlbabtlu 
OU M. destructor] D E] HL D. D] DE] U D] E] HL] RAPD E] UL UL GI. C] uutautuatutLutu Ll D DULL M. destructor | M. 
hordei l] B HH BU D HL D HO DL ELO HO DEBEO 

uguali 0000000 bU OO PCR-RFLPY RAPDU O O 

OU UU eees HOUUUDA HU D U D 00454-6296 2006[05-0822-07 



























































